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The new measurement of the anomalous magnetic moment of the muon by the Brookhaven AGS 
experiment 821 again shows a discrepancy with the Standard Model value. We investigate the 
consequences of these new data for neutralino dark matter, updating and extending our previous 
work [E. A. Baltz and P. Gondolo, Phys. Rev. Lett. 86, 5004 (2001)]. The measurement excludes the 
Standard Model value at 2.6a confidence. Taking the discrepancy as a sign of supersymmetry, we 
find that the lightest superpartner must be relatively light and it must have a relatively high elastic 
scattering cross section with nucleons, which brings it almost within reach of proposed direct dark 
matter searches. The SUSY signal from neutrino telescopes correlates fairly well with the elastic 
scattering cross section. The rate of cosmic ray antideuterons tends to be large in the allowed 
models, but the constraint has little effect on the rate of gamma ray lines. We stress that being 
more conservative may eliminate the discrepancy, but it does not eliminate the possibility of high 
astrophysical detection rates. 

PACS numbers: 95.35.+d, 14.80.Ly, 95.85. Pw, 95.85.Ry, 98.70.Rz 



> 

CO 

o 
o 

CN 
O 

^ ' 

6 ■ 
a: 

r/a . 



INTRODUCTION 



In early 2001, the Brookhaven AGS experiment 821 
measured the anomalous magnetic moment of the muon 
= (g — 2)/2 with three times higher accuracy than it 
was previously known [H . Their result disagreed with the 
Standard Model prediction at greater than 2.6a. How- 
ever, a sign error in the calculation of the hadronic light- 
by-light contribution to was discovered, reducing the 
discrepancy to 1.6a Recently, the same collaboration 
has released a result with much improved statistics ||, 
and there is again a discrepancy at the 2.6a level. Su- 
persymmetric particles can give significant corrections to 
a M |l| H H , thus the Brookhaven measurement is an im- 
portant constraint on supersymmetric models. There has 
been a substantial literature on this topic since the an- 
nouncement of the discrepancy j?], H , discussing various 
consequences of the older measurement. In this paper, we 
update the results of jjj concerning the implications of 
the Brookhaven data for supersymmetric cold dark mat- 
ter, assuming that supersymmetry is the only relevant 
physics outside of the Standard Model. 

There are two significant assumptions in our discus- 
sion. The first is that the Standard Model prediction for 
the muon anomalous magnetic moment is somewhat dis- 
puted, primarily in the hadronic contribution. This was 
clearly demonstrated in the sign error discovered in the 
last year. The hadronic error is a very significant part 
of the error budget when comparing the Brookhaven re- 
sults to the Standard Model. In fact it has been claimed 
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that the Standard Model errors have been significantly 
underestimated M , but this claim has been refuted [^0) . 
Furthermore, there are new evaluations of the hadronic 
vacuum polarization from firstly e + e~ — > hadrons indi- 
cating a larger discrepancy (3.6<r) and secondly hadronic 
tau lepton decays indicating a smaller discrepancy (1.3cr) 
|}| . The second caveat is that supersymmetry is only one 
of many possible scenarios providing corrections to a M at 
the weak scale. Theoretical prejudice tends to favor su- 
persymmetry, but other possibilities exist, summarized 
in Ref. 0. 



II. SUPERSYMMETRIC MODEL 

In the Minimal Supersymmetric Standard Model 
(MSSM) the lightest of the superpartners (LSP) is often 
the lightest neutralino. These four states are superposi- 
tions of the superpartners of the neutral gauge and Higgs 
bosons, 

X° = N n B + N 12 W 3 + N 13 H° + N 14 Hl (1) 

With R-parity conserved, this lightest superpartner is 
stable. For significant regions of the MSSM parameter 
space, the relic density of the stable neutralino is of the 
order fl x h 2 ~ 0.1, thus constituting an important (and 
perhaps exclusive) part of the cold dark matter (for a re- 
view see Ref. jllj ) . Note that Cl x is the neutralino density 
in units of the critical density and h is the present Hub- 
ble constant in units of 100 km s" 1 Mpc~ x . Large scale 
structure observations favor h = 0.7 ±0.1 and a matter 
density Qm = 0.3 ±0.1, of which baryons contribute a 
small amount n b h 2 = 0.02 ±0.001 @. CMB anisotropy 
measurements are consistent with this (summarized in 
e.g. PI), favoring Q M h 2 = 0.15 ± 0.05. We take the 
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TABLE I: The ranges of parameter values used in the MSSM 
scans of Refs. [h| |2o|. We use approximately 

80,000 models not excluded by accelerator constraints or the 
cosmological relic density bound (Q x h 2 < 0.25) before the 
measurement. Approximately 25,000 of these lie in the 



X < 0.25). 



unstable particles, all two-body tree-level annihilation 
channels of neutralinos, and coannihilation processes be- 
tween all neutralinos and charginos. We are in the pro- 
cess of including a complete treatment of sfermion coan- 
nihilations p3. 

Recent accelerator constraints are applied to each 
model in the database. Most important are the LEP 



bounds 1 25 on the lightest chargino mass (chargino and 
neutralino masses are tightly linked) 



88.4 GeV 

m + > < 
Xl 67.7 GeV 



|m x +-m x o|>3GeV ^ 
otherwise, 



range 0.05 < fl x h 2 < 0.25 as the cosmologically interest- 
ing region, basically a 2a constraint. Models where neu- 
tralinos are not the only component of dark matter are 
also allowed, so we separately consider arbitrarily small 
il x h 2 < 0.25. Even with very small relic densities, such 
models may be observable in astrophysical contexts [[l4j . 

We have explored a phenomenological variation of the 
MSSM with seven free parameters: the higgsino mass pa- 
rameter fj,, the gaugino mass parameter M%, the ratio of 
the Higgs vacuum expectation values tan (3, the mass of 
the CP-odd Higgs boson tua, the scalar mass parameter 
mo and the trilinear soft SUSY-breaking parameters Af, 
and A t for third generation squarks. All of our param- 
eters are fixed at the electroweak scale. Our framework 
is more general than the supergravity framework, in that 
we do not impose GUT unification of the scalar masses 
and trilinear couplings. In contrast to supergravity, this 
allows a highly pure higgsino LSP and its consequences, 
namely a SUSY spectrum that can be significantly more 
massive. For simplicity, we do apply the supergravity 
constraint on gaugino mass unification, though the re- 
laxation of this constraint would not significantly alter 
our results. As is typical, we assume R-parity conserva- 
tion, stabilizing the lightest superpartncr. (These models 
are described in more detail in Refs. |J, @ [f|.) 

To investigate the MSSM parameter space, we have 
used the database of MSSM models built in Refs. fl5|, 
101 [18, 



19, 20(1. Furthermore, for this work we have made 



special scans emphasizing large positive supersymmetric 
corrections to a^, Aa^(SUSY). The overall ranges of 
the seven MSSM parameters are given in Table [j| The 
database includes one-loop corrections for the neutralino 
and chargino masses as given in Ref. j2lj], and leading 
log two-loop radiative corrections for the Higgs boson 
masses as given in Ref. [ p2| . Supersymmetric contri- 
butions to the precision quantities a M and the b — > sj 
branching ratio are also included. The database con- 
tains the neutralino nuclcon cross sections and expected 
detection rates for a variety of neutralino dark matter 
searches. 

Crucial for studies of dark matter, the database in- 
cludes the cosmological relic density of neutralinos fl x h 2 , 
based on calculations in Refs. considering res- 

onant annihilations, threshold effects, finite widths of 



and on the lightest Higgs boson mass mh (which ranges 
from 91.5-112 GeV depending on tan (3) and the b — > sj 
branching ratio |^(| (DarkSUSY currently only imple- 
ments the leading-order calculation p7J). 



III. MUON ANOMALOUS MAGNETIC 
MOMENT 

Supersymmetric corrections to are surprisingly 
large, enhanced relative to typical weak-scale contribu- 
tions by the parameter tan/3 |6|. This fact makes 
these precision measurements enticing approaches for 
searching for supersymmetry. Typically, the supersym- 
metric corrections are given by 

A flM (SUSY)^14xlO- 10 ^^^) 2 tan/?, (3) 

where AfsuSY is the typical mass of superpartners. 

The new results of Brookhaven AGS experiment E821 
H for the anomalous magnetic moment of the muon, 
= (g — 2)/2, compared with the predicted Standard 
Model value are 

a M (exp) = 11 659 204(8) x 10~ 10 , (4) 
a M (SM) = 11 659 175(7) x 10" 10 , (5) 
Aa M = 29(11) x 10~ 10 . (6) 

This indicates a disagreement with the Standard Model 
at a 2.6cr confidence level. However, as mentioned in the 
introduction, there exist newer and conflicting evalua- 
tions of the Standard Model contribution which increase 
or decrease the discrepancy significantly. We will focus 
on the slightly older evaluation, but we will qualitatively 
discuss the effects of applying either of the new evalua- 
tions. To investigate the implications for the supersym- 
metric parameter space, we will assume that supersym- 
metry is the only source of corrections to a M outside of 
the Standard Model. Considering a 95% (2a) confidence 
region for the supersymmetric contribution, we accept 
the following range of Aa M (SUSY) 



7 x 10~ 1U < Aa^(SUSY) < 51 x 10" 
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(7) 
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We have used the full calculation in Ref. |jj to compute 
Aa M (SUSY) for the models in the database. 

The astrophysical phenomenology of neutralinos de- 
pends strongly on the ratio of gaugino and higgsino frac- 
tions, defined as 



Z n 
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(8) 



We plot this ratio against the neutralino mass for each 
model in the database in Fig. [l]. For clarity, the models 
have been binned along both axes. In the left panel, 
we only require that £l x h 2 < 0.25, and on the right we 
apply the more stringent constraint that the neutralino 
could make up all of the dark matter, VL x h 2 = 0.15 ± 0.1. 
Models allowed before the new Aa^(SUSY) constraint 
are plotted as crosses, and models respecting the new 
Aa^(SUSY) constraint are plotted as crossed circles. 

As has been discussed at length previously 0, ||, a 
Aa A1 (SUSY) bound that excludes zero from the positive 
side gives an upper limit on the mass of the neutralino, in 
this case 650 GcV. This is a large improvement over the 
cosmological bound based on the neutralino relic den- 
sity not being too large, an upper limit of 7 TeV p6| . 
However, if the Standard Model value is included in the 
allowed region by e.g. considering a 3cr confidence inter- 
val or a revised Standard Model calculation, there is no 
bound on neutralino mass. If the more favorable evalua- 
tion is considered, the mass bound becomes 450 GeV. 

If the neutralino has a large enough relic density to 
make up all of the cold dark matter, another effect ap- 
pears, namely that the neutralino can not be very purely 
higgsino-like in composition, requiring at least a 5% mix- 
ture (in quadrature) of gaugino states. Even without 
requiring neutralino dark matter, higgsino-like neutrali- 
nos are disfavored, with a maximum purity of 99.9% 

(Zg = 0.001). 

Orthogonal to the neutralino mass and composition 
but equally important to the value of Aa Al (SUSY) are 
the parameters tan/3 (the ratio of vacuum expectation 
values) and tuq (the scalar mass parameter). In Fig. ^ we 
plot these parameters for the database of models, again 
indicating the effects of the constraint on Aa M (SUSY). 
The constraint forces the scalar mass parameter to be 
small, but the upper bound increases with increasing 
tan/3. 



IV. ASTROPHYSICAL DARK MATTER 
SEARCHES 

There is a large community pursuing the goal of de- 
tecting dark matter particles, neutralinos especially, in 
various astrophysical contexts. The possibilities can be 
broken up along the lines of "direct" and "indirect" de- 
tection. Direct detection means detecting the rare scat- 
terings of neutralinos in our galactic halo with nuclei in 
a sensitive low background apparatus. Indirect detec- 
tion means detecting the products of rare annihilations of 



galactic neutralinos, such as antiprotons, antideuterons, 
positrons, gamma rays, and neutrinos. 

Perhaps the most promising of the astrophysical neu- 
tralino searches is the direct detection program. Ex- 
periments such as CDMS |§, DAMA §§], and EDEL- 
WEISS [3(J have pushed exclusion limits down to cross 
sections as small as 10~ 6 pb. As has been noted before, 
the neutralino-nucleon elastic scattering cross section ex- 
hibits a significant correlation with Aa A1 (SUSY) J?], |3l| , 
thus a large positive Aa /J (SUSY) is exciting for direct 
searches. Direct detection is promising even in the case 
where the neutralinos have a small relic density and thus 
are only a small component of the dark matter. In this 
case we perform a conservative rescaling of the galactic 
neutralino density as 



Px -> Px 



n x h 2 

0.25 



(9) 



where Sl x h 2 



0.25 is the current upper limit on the 
relic density. In the left panel of Fig. ||, we plot the 
spin-independent neutralino-proton scattering cross sec- 
tion, rescaled according to Eq. ^. The constraint due to 
Aa At (SUSY) is intriguing, as it bounds the rescaled cross 
section at around 10~ n pb. In the right panel of Fig. g, 
we perform no rescaling, and only consider models with 
cosmologically interesting relic densities. Here the mini- 
mum cross section is around 10~ 10 pb. The inlet at 100 
GeV and 10~ 9 pb is due to the lower limit Q, x h 2 > 0.05. 
These bounds indicate that there is considerable hope for 
the next generation of experiments, such as CDMS II and 
CRESST II §§. The latter bound is perhaps reachable 
by future experiments with one ton target masses such as 
GENIUS ||, CryoArray @, and XENON §|. Finally, 
it is important to note that in the case where the signifi- 
cance of the a M discrepancy is reduced, the lower bound 
on the cross section disappears. However, there is not an 
upper bound on the cross section. Large cross sections 
are still possible with Aa M (SUSY) consistent with zero. 

Neutrino telescopes such as at Lake Baikal [ 3q ] , Super- 
Kamiokande |37| , in the Mediterranean |3S ] , and the 
south pole [p9| are a promising technique for indirect 
detection. Neutralinos in the galactic halo scatter into 
orbits around the Earth or Sun, and can then rapidly 
sink to the cores of these bodies by additional scatter- 
ings, resulting in a large density enhancement. This can 
produce a detectable annihilation signal in neutrinos at 
high (GeV) energies. It is the capture rate that gov- 
erns the neutrino flux, and is strongly correlated with 
the neutralino-nucleon cross section. This places a lower 
bound on the detection rate, though at small neutralino 
mass there are threshold effects that remove it [fl9| . To 
illustrate, we plot the rate of neutrino-induced through- 
going muons from the Sun along with the unsubtractable 
background (from cosmic rays incident on the Sun's sur- 
face) in the top left panel of Fig. f| It is clear that the 
Aa M (SUSY) bound cuts away much of the undetectable 
parameter space, but not all of it. In the top right panel 
we repeat the calculation for neutrinos from the center of 
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FIG. 1: Gaugino/higgsino fraction versus mass for the lightest neutralino. In the left panel, we plot all models not excluded 
by cosmological arguments. In the right panel, only models with an interesting relic density are plotted. Crosses indicate 
models allowed before applying a constraint on Aa fl (SUSY), and crossed circles indicate models allowed after imposing the 
Aa M (SUSY) bound. 
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FIG. 2: Sfermion mass scale versus tan (3. As in Fig. gj, in the left panel, we plot all models not excluded by cosmological 
arguments and in the right panel we plot only models with an interesting relic density. Crosses indicate models allowed before 
applying a constraint on Aa M (SUSY), and crossed circles indicate models allowed after imposing the Aa M (SUSY) bound. The 
small "inlet" at tan (3 ~ 4 and mo ~ 400 GeV is due to the constraint on the Higgs boson mass. It is clear that the relic density 
cut has little effect on the allowed region. 



the Earth, where the prospects are much less promising. 

In addition, neutralinos can annihilate in the galactic 
halo. The relevant rates are quite small, but the enor- 
mous mass of the halo compensates, and the annihila- 
tion products may be detectable. Gamma rays propa- 
gate essentially freely, thus the expected rate is largest 
towards the galactic center where the dark matter density 
is largest. Charged particles are trapped by the galactic 



magnetic field and effectively diffuse, so these annihila- 
tion products would originate more nearby. 

The detection of the gamma ray lines from direct an- 
nihilations either to two photons, or to a photon and a 
Z boson would be a gold-plated signature of neutralinos 
in the galactic halo p8| . Gamma ray experiments such 
as the atmospheric Cerenkov telescopes (ACTs) VERI- 
TAS EfJ and MAGIC Ell, and the GLAST [El satellite 
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FIG. 3: Neutralino-nucleon elastic scattering cross section versus neutralino mass. In the left panel, we have only applied the 
upper constraint on relic density, and rescaled the effective cross section to account for a low galactic density of low relic density 
neutralinos. In the right panel we plot only those models that could account for all of the dark matter, and we do not perform 
a rescaling. The inlet at 100 GeV and 10 -9 pb is due to the lower limit Q x h 2 > 0.05. Crosses indicate models allowed before 
applying a constraint on Aa,j(SUSY), and crossed circles indicate models allowed after imposing the Aa,j(SUSY) bound. 



hope to detect these lines. Assuming that the galactic 
halo is an isothermal sphere with a 1 kpc core, we plot 
the reach of these experiments in the bottom left panel 
of Fig. ||. Note that with this assumption the emission 
enhancement from around the black hole at the galactic 
center is insignificant p3], so we neglect it. We notice 
that applying the Aa AI (SUSY) bound has little effect on 
the prospects for these experiments. It appears that the 
detection of the gamma ray lines is quite difficult. Other 
assumptions, including clumping of the dark matter or a 
lack of a central core lead to significantly higher predic- 
tions mi. 

The intensity of cosmic ray positrons p4| from neu- 
tralino annihilation is not much affected by the constraint 
on Aa M (SUSY). 

The final possibility we mention is that antideuterons 
may be an interesting annihilation product to search for 
p5[ . The background from mundane cosmic ray pro- 
cesses should be relatively smaller (a smaller fraction of 
the annihilation signal) at low energies than for antipro- 
tons. A signal may be detectable in experiments such as 
AMS_j4§] and GAPS J47|, as seen in the bottom right of 
Fig. [|. It is interesting that the Aa M (SUSY) constraint 
eliminates the models with the lowest rates, and further- 
more that the whole parameter space is covered for neu- 
tralino masses between 100 GeV and 500 GeV. Separat- 
ing a signal from the background with antideuterons may 
be difficult however. 



V. CONCLUSIONS 

In this paper we have discussed the recent confirma- 
tion of a discrepancy with the Standard Model of the 
anomalous magnetic moment of the muon ||, updat- 
ing and expanding the results of Ref. [Q. Assuming 
that supersymmetry is responsible for the discrepancy, 
we have investigated the consequences for astrophysical 
dark matter searches. We have confirmed that the con- 
straint significantly improves the prospects for direct de- 
tection experiments trying to measure the rare scatter- 
ings of galactic neutralinos. Neutrino telescopes are also 
helped by this result. The prospects for the detection 
of gamma ray lines from neutralino annihilations at the 
galactic center are not much affected. The prospects for 
detecting cosmic ray antideuterons as neutralino anni- 
hilation products are also significantly improved. In all 
cases, if the discrepancy disappears, there remain super- 
symmetric models with detectable rates for all of these 
experiments. 
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